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ABSTRACT 

I consider the equations of motion which follow from d'Alembert's principle for a general 
mechanical system in a space of N dimensions, constrained by a non-holonomic constraint 
which is linear and homogeneous in the generalised velocities. The variational equations of 
motion which follow for the same system by assuming the validity of a specific variational action 
principle, in which the non-holonomic constraint is implemented by means of the multiplication 
rule in the calculus of variations are also considered. It is shown that these two types of 
equations of motion are not compatible in a space of dimension iV > 3 if the constraint is 
genuinely non-holonomic. This means that these two types of equations of motion do not have 
coinciding general solutions. 
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1 Introduction 



Hamilton's principle for mechanical systems with non-holonomic constraints has recently been 
discussed by Flannery pQ. In particular a variational formulation of the equations of motion 
of a mechanical system was discussed both for holonomic and non-holonomic constraints. It 
was shown that while the equations of motion for a system with holonomic constraints can 
be obtained as variational equations, with the constraints being taken into account by the 
multiplication rule in the calculus of variations [2], the corresponding procedure with non- 
holonomic constraints leads to equations which differ from the equations of motion which follow 
from the well-known principle of d'Alembert. The questions discussed by Flannery are by no 
means new; they have been discussed in the literature in several papers. For a selection of 
references other than those given by Flannery, I refer to some of the references contained in a 
recent papei@ by Cronstrom and Raita [3J. 

Even though the equations of motion following from the principle of d'Alembert and from 
the variational action principle with non-holonomic constraints are different in form, there is still 
the possibility that the equations in question may have the same solutions. It was demonstrated 
by Pars [I], that this is not the case for a particular example with a specific non-holonomic 
constraint in a three-dimensional space. As such, this validates the assertion that the solutions 
to the two different types of equations of motion are different in general, at least in a space of 
three dimensions. It has only recently been proved, [3J, [5], that this is also the case generally 
in configuration spaces of dimension N > 3. Below I discuss an improved version of the proofs 
in [3J and [5], separately for N = 3 and N > 4. 

I consider a general autonomous system with a finite number of degrees of freedom, restricted 
only by reasonable smoothness conditions. For simplicity I consider only the case of one non- 
holonomic constraint, which is taken to be linear and homogeneous in the generalised velocities 
of the system. In principle the proof is valid when the configuration space of the system is a fairly 
general smooth iV-dimensional manifold, with N > 3. However, since the argumentation in the 
proofs of this paper for the most part is local, it is sufficient to use a co-ordinate formulation and 
consider only one appropriately chosen co-ordinate patch D q , with local co-ordinates designated 
as q = (q 1 , ...,q N ). 

In the case of two-dimensional systems there is nothing to prove, since every non-holonomic 
constraint of the kind considered in this paper can be reduced to an equivalent holonomic 
constraint when N = 2. 

The proof given in this paper implies definitely that the standard variational principle, 
which is valid for unconstrained systems, and which can be generalised to cover the case of 
holonomic constraints by using the multiplication rule in the calculus of variations, is not in 
general consistent with the principle of d'Alembert, if one uses the multiplication rule with a 
genuinely non-holonomic constraint. 

2 There are some regrettable misprints in this paper. The equations (17), (26) and (34) in Ref. [3] contain 
the redundant symbols = at the end of the equations. 
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2 The variational principles 



Consider an autonomous mechanical system with generalised co-ordinates q = (q 1 , ...,q N ), and 
velocities q = (q 1 , q N ). Throughout this paper I assume the the variables q = (q 1 ,...,q N ) 
belong to a chosen co-ordinate patch D q , which is a contractible domain in configuration space. 

The kinetic energy of the system is denoted by by T, and the generalised applied forces 
acting on the system are denoted by Qa,A — 1,...,N. It will further be assumed that the 
generalised forces can be expressed in terms of a generalised potential V(q, q), 

It is appropriate to introduce a Lagrange function L as follows, 

L(q,q) :=T(q,q)-V(q,q). (2) 

If the variables (q 1 , q ) are not constrained in any way, then, as is well known, the equations 
of motion of the system follow from the following variational principle, 

5S = 0, (3) 

where 

S:= JdtL(q,q). (4) 
The equations of motion are the following Euler-Lagrange equations, 

d ( dL(q,q)\ dL(q,q) 



dt \ dq A ) dq J 



0,A = 1,...,N. (5) 



I then consider a situation in which the variables (q 1 , q ) are constrained by a condition 
of the form 

*{q\...,q N ) = C, (6) 

where ^(q) is a smooth (continuously differentiable) function of its variables and C is a constant. 
The condition (jBJ) is a holonomic constraint. It is well-known that the equations of motion for 
the system also in this case can be obtained from the variational principle SS = 0, under the 
constraint ([6]). Using the multiplication rule in the calculus of variations [2] this leads to the 
following condition, 

5 J dt[L(q,q) + X^(q)] = 0, (7) 

which involves a Lagrange multiplier A. The Euler-Lagrange equations following from Eqn. (J7J) 
follows 

d ( dL(q,q) \ dL(q,q) 89(q) 

^{~W~)~~dq^~ = X -dq^~ i A = 1,...,N, (8) 
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The equations (jSJ) constitute a set of N second-order differential equations, which, under ap- 
propriate boundary conditions, and together with the constraint (jSj), are supposed to determine 
the quantities q 1 , q 2 , . . . , q N , as well as the Lagrange multiplier A. 

I then finally consider the case of a non-holonomic constraint, which is taken to be a linear 
and homogeneous in the generalised velocity variables q l , q N . The non-holonomic constraint 
is thus of the following form, 

N 

J2a A (q)q A = 1 (9) 
A=l 

where the N quantities a a (q), A = 1,2, . . . , N are given smooth (continuously differentiable) 
functions, which transform as covariant vector components under co-ordinate transformations. 

A metric g is involved in the notation in Eq. (j9]), which could also be written as follows, 

N 

£ 9AB(q)a A (q)q B = 0. (10) 
A,B=1 

Tensor- indices are thus lowered or raised with the metric tensor gAB^q) and g AB (q), respectively, 
where 

| ito(5 ) S -( 9 ) = ff . C :={j l A A = l (ID 

for A,C =1,2,... ,N. 

It is clear that that not every component of (^(g) can vanish identically in D q , since other- 
wise there would not be any non-holonomic condition at all in the problem. 

The question is then whether one can use the analogue of the equations (J7|) for a non- 
holonomic constraint of the form ([9]), simply by enforcing the constraint (Q by means of the 
multiplication rule throughout a domain which is supposed to contain all the competing paths 
in the variational principle. 

I now prefer to use the special notation 

L (q,q) = T(q,q)-V(q,q). (12) 

This special notation is used as a reminder of the fact that the quantity L (q, q) that one should 
be aware of the potential danger involved in using the quantity L Q as a Lagrangian in Hamilton's 
principle when the constraints are genuinely non-holonomic. In the absence of a constraint of 
the form (j5J), the quantity L (q, q) would naturally be the Lagrangian of the system. 

The following condition is a straightforward analogue of the variational condition (J7J) , 



5 / dt 



N 



L (q,q) a A (q)q J 

A=l 



0, (13) 



which involves a Lagrange multiplier The variational principle underlying the condition 
(fl3"j) is, as such, a possible variational principle for certain types of problems, but it does not 
necessarily give rise to appropriate equations of motion in Newtonian classical mechanics. 
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The condition (1131) gives rise to the following Euler-Lagrange equations, 
d (dL {q,q)\ dL Q {q,q) . * 

it ) - = * aA{q) + " £ ABiq)q ! = 1; ■ ■ ■ ' N > (14) 

where the quantities Mab(?)> A, B = 1.2. ... ,N are given as follows, 

Mab{q) '■— ~~QyB dq~A~ J AS = l,...,iV. (15) 

The variational equations of motion (|14p are precisely those which have been proposed from 
time to time in the literature, as appropriate variational equations of motion for mechanical 
systems with non-holonomic constraints of the type considered in this paper. 

It should be noted that the equations (IT4"|) and the constraint admit a first integral for 
the system, 

E -=2^1 —HTa L (q,q), (16) 

A=l u 1 

where E is the constant energy of the system. 

Consider now the special case in which the constraint (Q is integrable. This means that 
there exists a function ^f(q), say, such that 

a A (q) = ^(q),A=l,2,...,N. (17) 

Then the constraint (|9]) is in fact a holonomic constraint of the kind given in Eq. (J6]). 

Inserting the expressions (1171) in the equations f[T4l one obtains the following equations, 



d ( dL(q,q) \ dL(q, q) d9(q) 

Jt\-^)-^q-^ = ^^' A=1,...,N, (If 
which naturally agree with the equations OH]) when one makes the change of notation fi — y A. 



3 Genuinely non-holonomic constraints 

It was noted above that if there exists a function ^(q) such that the condition ( ITT)) holds true, 
then the constraint (Q is equivalent to a holonomic constraint of the kind given in Eq. ()6]). If 
the condition (fTTj) holds true, then necessarily, 

M AB (q) = 0, A, B = 1, 2, . . . , N, q E D q , (19) 

where the quantities Mab,A,B = 1,2,...,N were defined in Eqns. ( fTST) . Conversely, in a 
contractible domain, the integrability conditions ( fT9l) are also sufficient for the integrability 
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of the constraint ([9]). Thus, the constraint fl9]) is integrable, and therefore equivalent to a 
holonomic constraint, if and only if the integrability conditions f fl9|) are fulfilled. 

In the case at hand, the constraint is equivalent to the following constraint, 

N 

$(<?) £ a A (q)q A = 0, (20) 

A=l 

where <&(q) is some appropriate smooth function, which does not vanish identically. It may 
be that one can choose the function <&(q) in ( 1201) in such a manner that this constraint is 
integrable. If that is the case one says that $(g) is an integrating factor. Replacing a A by 
$>a A ,A = 1,2, ...,N in the integrability conditions (fl9|) . one finds the following necessary 
conditions for the existence of an integrating factor, 

a A (q)M BC {q) + a B {q)M CA {q) + a C {q)M AB {q) = 0, A, B, C = 1, 2, . . . , N. (21) 

The number of independent integrability conditions of the type (|2ip is N c , where [8], 

N c = ±(N-l)(N-2). (22) 

In three-dimensional space (N = 3) there is thus only one integrability condition of the type 
(j2ip . namely the following condition, 

ai(q)M 23 (q) + a 2 (q)M 3l (q) + a 3 (q)M 12 (q) = 0. (23) 

In the case of a contractible domain the integrability conditions (121 j) are also sufficient for the 
existence of an integrating factor <&(q). 

Thus, if there exists an integrating factor $(<?), then the original constraint (|9l) , which has 
the appearance of a non-holonomic constraint, can be replaced by an equivalent holonomic 
constraint by making use of the integrating factor. 

In what follows it is important to consider only such constraints which are genuinely non- 
holonomic. A genuinely non-holonomic constraint, linear and homogeneous in the velocity 
components q A , A = 1, 2, . . . , A, is a constraint of the form Q which is neither integrable as 
such, nor integrable by means of an integrating factor. Thus, for a genuinely non-holonomic 
constraint neither the integrability conditions ffl9|) nor the integrability conditions (121]) are 
fulfilled. 



4 The d'Alembert equations with constraints 

The principle of d'Alembert (see e.g. the classical texts by Goldstein |6] or Whittaker [7]) gives 
the following equation, 



N 

i x \dt \dq A J dq 1 
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where the quantities 5q A are virtual displacements of the system. The meaning of the symbols 
T and Qa, respectively, is the same as before, i.e. T is the kinetic energy and Qa is the A:th 
component of the generalised external force imposed on the system. If the virtual displacements 
Sq A ,A = 1,...,N, are independent, then Eq. (|24|) results in the ordinary d'Alembertian 
equations of motion, 

Consider the generalisation of the discussion above to systems with a constraint of the form 
([9]) given above. The derivation given below of the equations of motion for this system, with 
a potentially genuinely non-holonomic constraint, can be found in the textbook by Whittaker 
0. 

Implement the constraint fl9]) by regarding the system to be acted on by both the external ap- 
plied forces Qa, A = 1, . . . , N, and by certain additional forces of constraint Q' A , A = 1, . . . , N, 
which force the system to satisfy the condition (Q. Equation (pHl) is then replaced by the 
following equation 

d_(dT_\ dT 
Jt ~ dq 

In Eq. (1261) the virtual displacements 8q , A = 1, . . . , N, can now be regarded as independent 
Thus one obtains the equations of motion, 



N 



E \tA^a)-^a-Q a -Q'a\Sc1 A = 0, (26) 

A=l 



it(w)-w^ QA+Q '- A ^'-' N - (27) 

The forces of constraint, Q' A , A = 1, . . . , N, are a priori unknown, but they are such that, in 
any instantaneous displacement 5q A , A = 1, . . . , N, consistent with the constraint ([9]), they do 
no work. The non-holonomic constraint (Q implies the following condition on the possible 
instantaneous displacements Sq A , A = 1, . . . , N, 

N 

53 a A (q)5q A = 0. (28) 

A=l 

For any instantaneous displacements 5q A , A = 1, . . . , N, which satisfy the condition (128]) . the 
work 5W done by the constraint forces Q' A , A = 1, . . . , N must be equal to zero, i.e., 

N 

SW' := Q'a^ A = 0. (29) 

A=l 

The conditions (|28|) and (|29|) together imply that 

Q' A = Xa A (q),A=l,...,N, (30) 



where the quantity A is a time- dependent parameter to be determined. Eqns. ( 1271) have thus 
been reduced to the following form, 
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These N equations of motion are consequences of the principle of d'Alembert. They will be 
called d'Alembertian equations. One should still add the equation of constraint (jUJ) to the 
equations of motion above. There are thus altogether N + 1 equations for the determination of 
N + 1 quantities q (t), A — 1, . . . , N, and A(t), when appropriate boundary conditions for the 
quantities q 1 , q N and q 1 , q N are given. 

One should observe that in the argument above, it is not required that the constraint 
equation be in force under general variations qj —> qj + 5qf, the constraint (J2J) is only 
imposed on the actual motion of the system. 

Assuming that the generalised force components Qa, A = 1,2, . . . , N can be expressed in 
terms of a generalised potential V(q,q), as in Eq. (CE]), one can rewrite the d'Alembertian 
equations ( 13TT) as follows, 

d (dL (q,q)\ dL (q,q) 

where the quantity L (q, q) has been defined in Eq. ( Tl2l) . 

The equations fl32|) will also be referred to as d'Alembertian equations in what follows. 

The d'Alembertian equations of motion ( 132]) admit a first integral, namely the expression 
( TT6]) . This is the same first integral as the one obtained from the variational equations ([1] 



Non-equivalence of the principle of d'Alembert and 
the variational action principle with genuinely non- 
holonomic constraints 



I now consider the d'Alembertian equations of motion fl32|) and the variational equations (JT^ 
respectively, as initial value problems. This means that the co-ordinate patch in which the 
quantities q l ,q 2 , . . . ,q N are local co-ordinates, is assumed to be chosen so that it includes the 
initial value point go = (<?d, Qoi ■ ■ ■ > ^oO- 

The question is whether the equations f )32|) and (j!4p . respectively, can have the same solu- 
tions for . . . ,q N {t) in general, despite the fact that these equations are not identical. It 
will be shown below that this is not the case. 

Naturally, the functions L Q (q,q) and aA^q) entering into the equations fl32|) and ffT^j) . re- 
spectively, will have to satisfy appropriate smoothness conditions in order that these equations 
may have solutions. I will not enter into a discussion of such smoothness conditions, but rather 
assume that the equations ( 1321) and respectively, have e.g. C 2 -solutions in some appro- 
priate time-interval, for given initial values for the co-ordinates q A in the domain D q and for 
given initial values for the velocities q A , A = 1, . . . , N, 

~q A (t)] t=to = q A , [q A (t)] t=to = q A , A = 1, . . . , N. (33) 
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The initial values q$ and q^ at t = to are free parameters within an appropriate region of the 
configuration- and velocity space, except for the restriction 

N 

°a(<?o) Qo = 0- (34) 

A=l 

The condition (134]) is a consequence of the non-holonomic constraint (Q. 

After these considerations concerning the initial values I will prove a theorem on the incom- 
patibility of the d'Alembertian and variational equations of motion, respectively. 

5.1 The incompatibility theorem 

Consider the d'Alembertian and variational equations of motion ( 132]) and (}T4"j) . respectively, 
which involve the genuinely non-holonomic constraint (Q. If N = 3 the d'Alembertian and 
variational equations of motion are not compatible in the sense that they do not have coinciding 
solutions. If iV > 4 the d'Alembertian and variational equations of motion are not compatible 
in the sense that they do not have coinciding solutions with arbitrary general initial values ( )33l) . 
which satisfy the condition (|34|) . 

It is assumed that the d'Alembertian equations (I3"2"|) have unique smooth (e.g. C 2 ) solu- 
tions (q 1 (t) , . . . , q (t)) in an appropriate time-interval, satisfying a genuinely non-holonomic 
constraint of the type (]9]), and the general initial value conditions ( 133]) and (134"]) . I call such 
solutions general solutions in what follows. 

The method of proof is by reductio ad absurdum, i.e., one makes the assumption that 
the variational equations (]T4"]) have solutions which coincide with the general solutions of the 
d'Alembertian equations ( 132|) . and shows that this assumption leads to contradictions. 

The proof given below is designed in such a manner that it is independent of any specific 
properties of the kinetic energy T and generalised potential V, respectively, for the problem 
under consideration. 

5.2 Preliminary considerations 

Assume now that the equations (|14p and f ]3"2"]) have coincident solutions. By subtracting Eqns. 
(fi4|) from Eqns. ( 132]) . one obtains the following equations, 

TV 

(A - fi) a A (q) = fM J2 M Mq) q B = 0, A = 1, . . . , N, (35) 

B=l 

which have to be satisfied by the general solutions (q l , . . . , q ) of the d'Alembertian equations 
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It should be noted that one must necessarily have \i ^ in the equations f|35|) above, 
since otherwise the variational equations (fT4l would contain no reference whatsoever to the 
constraints ([9]). Further, if fj, = 0, it follows from the equations ( 135]) that also the following 
conditions must hold true, 

Xa A eeO, A=1....,7V. (36) 
The conditions (|36l) above imply that either 

A = 0, (37) 

or 

a A (q) = 0, A=l....,iV. (38) 

If the conditions f )38|) are true, there are no constraints to be considered so the the basic 
question under consideration disappears. Finally, if the condition ( 1371) holds true in addition 
to the condition /i = 0, then the d'Alembertian equations of motion ( 132]) and the variational 
equations of motion (fl4j) are trivially identical, so the the basic question under consideration 
disappears also in that case. Hence one must have /i ^ 0. 

Using the notation 

T := (39) 
li 

one can write the conditions ( 1351) in an equivalent form as follows, 

N 

Y,M AB (q)q B = Ta A (q), A = 1,...,N, (40) 

B=l 

The tensor equations (140 p constitute a set of N linear algebraic equations in the variables 
q A , A — 1, . . . , N, for any fixed value of q := (q 1 , q 2 , . . . , g w ) G D q . It should be noted that the 
quantities M A B(q) in Eqns. ( 140]) are anti-symmetric upon an interchange A -h- 5 of the indices 
A and fi, 

M AB (g) = -M BA (g), A B = 1, 2, . . . , N. (41) 
Thus the matrix-like quantity M(q) := (M A B{q)) is skew. 



5.3 The case N = 3 

I consider the case N = 3 separately, since in that case the analysis of the tensor equations (14T)]) 
is particularly simple. Thus, for N = 3 the equations (1401 read as follows, 

3 

Y,M AB (q)q B = Ta A (q), A =1,2,3. (42) 

5=1 

It should be noted that since N = 3 and M (q) is skew, one necessarily has 

det M(q) = 0. (43) 
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The quantity T in Eqns. (142]) is a so far unknown parameter, but one must have either 
T = or r ^ 0, 

Consider first the case T = 0. Then the velocity components (g 1 , q 2 , q 3 ) satisfy the homo- 
geneous equations 

^M AB (q)h B (q)=0, A = 1,2,3. (44) 

B=l 

The homogeneous equations (144 p have the following non-trivial solutions, 

(h\q),h 2 {q),h 3 (q)) := a{M 23 {q) , M 31 (q) , M X2 {q)) , a^O. (45) 
However, the solutions fj45|) for (g 1 , q 2 , q 3 ) must also satisfy the constraint ([9]) with N = 3, i.e., 

N 

£a A (g)g A = 0. (46) 

A=l 

This constraint then implies that 

a 1 (q)M 23 (q) + a 2 (q)M 31 (q) + a 3 (q)M 12 (q) = 0. (47) 

The equation ( l4Tj) is nothing but the condition ( 123|) . which implies that the constraint ( 146]) is 
integrable by means of an integrating factor. This is a contradiction, since by assumption the 
constraint ( 146]) is genuinely non-holonomic. 

If T ^ 0, the equations ( 1421) are inhomogeneous equations. In order that the inhomogeneous 
equations (142]) may have a solution it is necessary that the quantity a = (ai, a 2 .a 3 ) on its right 
hand side be orthogonal to any non-trivial solution (145]) of the homogeneous equation ( 144]) . 
This also leads to the condition (147]) as a necessary condition for the solvability of Eqns. ( 142]) . 
This implies the same contradiction as in the homogeneous case. 

It has thus been proved that the d'Alembertian equations of motion ( 132]) and the variational 
equations of motion (1T4]) in the case iV = 3 can not have coincident solutions if the constraint 
(146 p is genuinely non-holonomic. 



5.4 The cases N > 4 

I will now consider the linear equations (140]) for any finite dimension iV > 4, at some fixed 
point q = (q 1 , q 2 , . . . , q N ) G D q . 

The equations ( 140]) . which are necessary consequences of the assumption that d'Alembertian 
equations of motion ( 132]) and the variational equations of motion ( 1T4]) have coincident solutions, 
are valid at any given fixed point q G D q . 

It is now appropriate to analyse the equations ( 140]) in more detail. Contract Eqns (140]) with 
the quantity (a 1 (g), a 2 (q), . . . , a N (q)). This leads to the result 

N 

£ a A { q )M AB {q) q B = T\\a{ q )\\\ (48) 

A,B=1 
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where the norm ||a(<7)|| is defined by the expression 



N 

l«(g)l| 2 =E« A (^W(?). (49) 

A=l 



Let R(q, q) be defined as follows, 



N 

R(q,q):= E a A (q)M AB (q) q B . (50) 
A,B=l 



Trading the quantity Y for R(q, q) in Eqns. (HO"]) one obtains, 

Y. M A D (q)q D = " A !f} R }lS = 1,2, . . . , AT, 



_D = 1 



K<z)IP 



(51) 



where the free index A is raised for convenience. Contracting Eq. (l5Tj) with the quantity 
J2b=i MabQ b , and using the anti-symmetry condition fT4"Tj) . one gets finally the following ex- 
pression for the quantity R(q, q) defined in Eq, floTJ]) . 

N 

R 2 (q,q) = -\\a(q)\\ 2 J2 (M 2 ) BD (q)q B q D , (52) 

B,D=l 

where 

A 7 " 

(M 2 ) BD (g):=^M Byl (g)M A D (g), fl,£> = 1,2, . . . ,JV. (53) 

A=l 



From the anti-symmetry condition f HTj) follows also that the quantity M 2 is symmetric in its 
lower indices, 

{M 2 ) BD {q) = (M 2 ) DB (q), B, D = 1, 2, . . . , N. (54) 

Incidentally, Eq. f )52|) shows that the following quadratic form Q(2){x) in real variables x A ,A = 
1, 2, . . . , N, is negative semi-definite, 

N N 

Q{2){x)-= £ E M B A(q)M A D (q)x B x D < 0. (55) 

B,D=1 A=l 

This is made even more explicit by a reshuffling of the summation indices in Eq. ( 15 5p . 

N / N \ f N \ 

Q(2)(x) = - E /°(?) E ^ B (g)x B E ^CD)g)^ D = -\\M(q)x\\ 2 . (56) 

A,C=1 \B=1 / \D=1 / 

Thus the quadratic form Q(2)(x) vanishes if and only if 

N 

E M AB (q)x B = 0. (57) 
The results given above will be used subsequently. 
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Applying the formula fl56|) to the expression fl52|) one gets finally the following useful ex- 
pression, 

R(q,q) = + \\a(q)\\\\M(q)q\\. (58) 

There is an ambiguity involving a sign ±1 in arriving at the expression ( 158]) along the route 
indicated above. However this sign-ambiguity is of no consequence for the remaining consider- 
ations in this paper, so for simplicity I shall stick to the sign- factor +1 in the expression (|58l 
in what follows. 

Inserting the expression (1581) for R(q, q) in Eqns. ( l5~Tj) one finds the following equations, 

E MAB{q)q B = "if^n 11 a A (q),A = 1,2,...,N. (59) 

The tensor equations (159|) are more precise versions of the equations ( 140|) . which constituted 
the starting point in the proof of the non-compatibility theorem. 

It is necessary to consider the detailed structure of the solutions q A ,A = 1,2, ...,N, of 
Eqns. (159|) in order to proceed with the proof of the non-compatibility theorem enunciated 
in Subsection 15.11 For this purpose one needs certain basic facts concerning skew matrix-like 
quantities such as the anti-symmetric tensor fields (Mab((/)) in the equations above. 



5.5 Digression on skew matrices (Mas) 

In this subsection I will discuss such properties related to the skew matrix-like quantities 
M A B{q), A, B = 1,2, . . . , N, encountered previously, which are needed for finding solutions 
of Eqns. pD for q A ,A = 1, 2, . . . , N. The components g AB {q), A, B = 1, 2, . . . , N of the 
metric tensor g enter in an essential way in this discussion. 

The co-ordinates q enter as parameters in the discussion below, which deals with exclusively 
algebraic properties. In this discussion the co-ordinates q may either considered to correspond 
to a fixed point q G D q , or else to be confined to an appropriate reducible sub-domain D® C D q , 
which contains the given fixed point go- Eventually the fixed point go will be identified with 
an initial value in co-ordinate space for the d'Alembertian and variational equations of motion 
under consideration in this paper. 

Consider the following quadratic form Gm)(x) in certain real variables x A ,A = 1, 2, . . . , N, 

N 

G {2) {x) := E 9AB(q)x A x B . (60) 

,4,5=1 

It follows from the assumed properties of the metric g(q), that the quadratic form G( 2 )(x) is a 
real positive definite quadratic form. 

I also recall the negative semi-definite quadratic form Q(2){x) defined above in Eq. ( )55|) . 
involving the quantities (M 2 ) BD = (M 2 ) DBl B, D = 1,2, . . . , N defined in Eq. 
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It is known that one can effect a simultaneous reduction of two quadratic forms, of which 
one is a positive definite form such as the form (1601) above, and the other is merely a real 
quadratic form, such as the form ( 155]) above, to sums of perfect squares [9|. This leads to the 
following eigen- value problem for eigen- vectors e A (q), A = 1, 2, . . . , N, 

N 

E {(M 2 )AB(q) ~ X(q) 9AB(q)) e B (q) = 0. (61) 

B=l 

The equation (1611) has non-trivial solutions only if the parameter A satisfies the following eigen- 
value equation, which is similar to the characteristic equation for eigen-values in ordinary linear 
algebra, 

det ((M 2 ) AB (q) - \{q) g AB {q)) = 0. (62) 

Under the conditions described above for the quadratic forms G( 2 )(x) given in ( 16"U|) . and the 
quadratic form Q(2){ x ) given in ( 153]) . the equation (|6"2"|) has N eigen-value solutions \ u (q), v = 
1,2, . . . , N. The corresponding eigen- vectors e u B (q) satisfy the following equations, 

N 

E {(M 2 ) AB (q) - \ u (q) g AB (q)) e B (q) = 0, u = l,2,...,N. (63) 

B=l 



It is convenient to consider the following variant of Eqns. (163]) . 

N 

E ((M 2 ) A B (q) - \ u (q) g A B (q)) e vB {q) = 0, v = 1, 2, . . . , N. (64) 

B=l 

The equations (|63]) and (164]) are perfectly equivalent; the equations (1641 are obtained by simul- 
taneously raising and lowering the summation indices B inside the equations (163]) . which is a 
straightforward and legitimate operation. 

From Eqns. (16"4"]) follows that the eigen-values A„(g), v = 1, 2, . . . , N, also satisfy the follow- 
ing variant of Eq. (152]) , 

det((M 2 ) A B (q)-\ u (q)9A B (q))=Q- (65) 

Here the enumeration of rows and columns in the determinant is as follows: the lower indices 
A = 1,2, ... ,N enumerate the rows, and the upper indices B = 1,2, ... ,N enumerate the 
columns. 

The equation (165]) has an advantage over Eq. (]62]) . namely that the metric g is only 
involved in the expressions (M 2 ) A B (q), A, B = 1,2, . . . , N, in Eq. (]63]) . This is due to the 
fact that the quantities g A B (q), A, B = 1,2, ... ,N always have the fixed numerical values or 
1, independently of the co-ordinates q, 

b / \ J 1 if A = B 
9A (?):=| o iiA _, Bj (66) 

for A,B = 1,2, ...,N. 
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It should also be noticed that it follows readily from Eq. fl65|) that the eigen- values X u (q), v = 
1, 2, . . . , N, are scalars, i.e. invariant under co-ordinate transformations. 

I finally rewrite Eqns. ( I64p as follows, 

N 

£ (M 2 ) A B (q) e uB (q) = X u (q) e uA (q), v = 1, 2, . . . , N, A = 1, 2, . . . , N. (67) 

B=l 

It should be noted that the eigen-vectors e v A{q), A = 1,2, . . . , N are orthogonal vectors, as 
will be shown below. 

If the eigen- values \ u (q), is = 1,2, N are non-degenerate, this follows immediately from 
a consideration the equation (l6Tj) for two different indices n and v, say. Contract the equations 
(!67l) with e^g). Interchange ft and z/ in the equations obtained in this way, and subtract 
one set of equations from the other set of equations in the pairs of equations obtained by the 
contraction procedure. Then one obtains, 

(A M (g) - A„(g))(e M (g), e vB {q)) = 0, (68) 

where 

TV 

Mq),e v {q)):= £ g AE \q)e, A {q)e uB {q). (69) 

A,B=1 

The abbreviated notation on the left hand side of the expression (16"9~1) for the inner product of 
vector-quantities such as the eigen-vectors e^g), A = 1, 2, . . . , N, will be frequently used in 
what follows. 

If the eigen- values are non-degenerate, i.e., if A M (g) ^ X u (q) when ^ v, then the or- 
thogonality of the corresponding eigen-vectors e^(q) and e u (q) follows from the equations (!68l) . 
Normalising the eigen-vectors to 1 in the inner product (|69|) . one then finally obtains the fol- 
lowing ortho-normalisation conditions, 

for fj,,u = 1,2,..., N. 

If there are degenerate eigen- values, the ortho-normalisation conditions are imposed as a part 
of the definition of the corresponding eigen-vectors. Thus the ortho-normalisation conditions 
([70]) are valid in all cases. 

Using the ortho- normalisation conditions flTOj) one readily obtains the eigen- values X u {q), v = 
1,2,..., N, from Eqns. (1671). 

N 

K(q)= E e u A (M 2 ) A B (q)e l/B (q))^-\\M(q)e u (q)\\ 2 ^ = ^---,N. (71) 

.4,73=1 

In arriving at the final expressions above for the eigen- values \ v ,v = 1,2, . . . ,N , use has been 
made of the result given in Eq. (1561) . 
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From the expressions (ITT]) follows that all the eigen- values A M ,/i = 1,2,...,N are non- 
positive, 

A„(g) := -^(g) 2 , (72) 

where the real quantities K u (q), v — 1, 2, . . . , , N, can be chosen to be non-negative without loss 
of generality. Comparing Eqns. (J7T1) and ( 1721) one then obtains, 

^(g) = ||M(gK(g)||, */ = 1, 2, . . . , N. (73) 

From Eq. (J73~l) one gets an upper bound on the quantities k„, v — 1, 2, . . . , N, in terms of the 
norm ||M(g)|| of the matrix M(q), 

\\M(q)\\ = sup M=1 \\M{q)x\\. (74) 

The quantities k u ,u = 1,2, . . . , N, thus satisfy the following inequalities, 

<«„(</)< | |M(?)| |, v=l,2,...,N. (75) 



From Eq. ( !7T|) follows that an eigen-value \ v (q) equals zero if and only if 

N 

J2M AB (q)e u B (q) = 0, A = 1, 2, . . . , N. (76) 

B=l 

The number of non- vanishing eigen- values is thus determined by the rank of the matrix M(q). 
Since M(q) is skew, the rank of M(q) is necessarily an even integer, 2p, say, where p > 1, since 
M(q)^0. 

The eigen- values will now be enumerated in such a manner that the last N — 2p eigen- values 
X u ,v — 2p + 1, . . . , N are equal to zero. Thus, 

JV 

E^Afl(?)e/(?) = 0, A = l 1 2,...,N;is = 2p + l 1 ...,N. (77) 

B=l 

The subspace ker(M(q)) is thus spanned by the eigen-vectors e u (q), v = 2p + 1, . . . , N, the 
components of which have been labelled by upper or indices A = 1,2, . . . , N, in Eqns. (1771) . 
One may also use lower indices A, B, . . ., whenever it is convenient, by using the standard rules 
for lowering or raising such indices with the metric tensor gAB^q) or g AB (q), respectively, where 
A, B — 1, 2, . . . , N. Corresponding to the enumeration in Eqns. (ITTj) one recognises that 

K v (q) = 0, u = 2p+l,...,N. (78) 

If M(q) is regular there are no non-trivial solutions of the equations (ITTj) . This can happen 
only if iV is even and when 2p = N. 

Generally, for v = 1, . . . , 2p, the quantities K u (q) are positive, 

Ky(q) > 0, v= l,2,...,2p. (79) 
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A remark on the classification of a skew matrix M(q) according to its rank may be in order. 
Since M(q) depends on q, also the rank 2p is in principle dependent on q. If considerations 
related to a skew matrix M(q) involving only one prescribed fixed point qo, say, no problems can 
occur by the omission of the argument q in the quantity p. This notation is used here for the 
sake of simplicity. Strictly speaking, however, one should use the notation p(q) in stead of only 
p. The rank 2p(q) takes values in the set of even positive integers, and therefore the rank 2p(q) 
ought not to change from its value at a prescribed point qo if one considers neighbouring points 
q sufficiently near the prescribed point qo- This also justifies the omission of the argument 
q from the symbol 2p, which then can be considered to be a valid characterisation of M(q) 
at least in an appropriate sufficiently small domain containing the prescribed point qo as an 
interior point. 

I will now construct a new basis in the JV- dimensional linear space spanned by the eigen- 
vectors e u (q), v = 1, 2, . . . , N. This basis is used to solve the important equations (15Uj) . The 
construction parallels closely a similar construction made by Greub [TO] for skew mappings in 
ordinary Euclidean space. 

Define vectors & 2 i/-i(<z)> ^ = 1, • • • as follows, 

K-Aq) ■= e„ A (q), v = i> • • • ,p ; A = 1, 2, . . . , JV. (80) 

Similarly, let vectors b 2u {q), v — 1, . . . ,p, be defined by the equations, 

N 

b 2 uA(q) ■= ^u(q)- 1 E MAB(q)K-i B (Q), " = h • • • >P I A = h 2 , . . . , N. (81) 

B=l 

In view of the conditions ( I79p the definition is a possible definition. 

I rewrite the eigen-value equations f lBT|) with —n u (q) 2 in stead of \ u (q), thus 

N 

E ( M2 )A B (q) euB{q) = -kMY (?) e„ A (q), v = 1, 2, . . . , N, A = 1, 2, . . . , JV. (82) 

B=l 

It is now a simple matter to show that 

N 

E (M 2 ) A B (q) b^q) = -K v {qf b^ A (q), ^ = 1, 2, . . . , 2p, A = 1, 2, . . . , JV. (83) 

B=l 

When 2p < N one also defines, 

b vA {q) ■= e uA (q), v = 2p + 1, . . . , N, , A = 1, 2, . . . , JV. (84) 

Then, 

N 

J2(M 2 ) A B (q)K B (q)=0, v = 2p + l,...,N;A=l,2,...,N. (85) 

B=l 
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From the definition flHTj) follows that 

N 

X; M AB ( g )^_ 1 B ( ? ) = ^(g)6 2 ^( g ), )I / = l,...,p : A=1,2,...,N. (86) 
Using the equations ( 1851) it is also simple to show that 

N 

X M AB (g)^ B (g) = -^( g )6 2l ,_. 1A (g), , 1/ = 1, . . . ,p : A = 1, 2, . . . , N. (87) 

B=l 

It is finally stated without proof, that the basis vectors b^(q),fi = 1,2, ... ,N, defined above, 
are orthonormal in the inner product (169|) . 

%(q),K(q)) = S^, p,v=l,2,...,N. (88) 

The proof of the ortho-normalisation conditions ( 1881) is fairly simple, and is left to the interested 
reader. 

The vectors , p = 1,2, . . . , N thus also form an orthonormal basis in the space spanned 
by the eigen- vectors e M , [i = 1, 2, . . . , N of the matrix M 2 (q). 

Incidentally, the results given above imply that there are only p non-zero eigen-values 
^v(q) — —K u (q) 2 related to the matrix M 2 (q) in stead of 2p such values. This means that 
the non-zero eigen-values of M 2 (q) are at least two- fold degenerate. 



5.6 The solutions q 

I now return to the important equations ( 159|) . These equations are valid for q e D q , but will 
later only be considered at a fixed point qo, which is identified as the initial value point for 
the d'Alembertian and variational equations of motion f)32p and (Il4p . respectively. The initial 
velocity at go will be denoted by qo. 

For simplicity of notation, I will for the time being consider the equations at an arbitrary 
fixed point q e D q . 

Contracting the equations ( 15 9 j) with an arbitrary basis- vector b u (q), with v in the range 
2p + 1, . . . , N, one obtains, 

(K(q),a(q)) \\M(q)q\ | = 0, v = 2p + 1, . . . , N. (89) 

From Eq. (|89|) follows that either, 

||M(g)g||=0, (90) 

or, 

\\M(q)q\\^0, (91) 
in which case the following conditions must also be fulfilled, 

(K(q),a(q))=0, v = 2p + 1, . . . , N. (92) 
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5.6.1 The case ||M(g)g|| = 0. 

Consider first the case (19U1) with 2p < N. Thus g G ker(M(q), or, equivalently, 

N 

Q A = E l»K A {q), A = 1,2,..., N. (93) 

i/=2p+l 

The parameters j v , v = 2p + 1, . . . , N, in the expressions fl95|) are restricted by the non- 
holonomic condition (Q. Thus, 

N 

= (a(g),g) = £ 7„(a(q), b„(q)). (94) 

v=2p+l 

I now consider the results above at the fixed initial value point go- 
The equation (1931 then reads as follows. 

N 

4o A = E 7o*&/(gb), A =1,2,..., iV. (95) 

i^=2p+l 

The parameters 7o„, f = 2p + 1, . . . , , N, are restricted by the condition (|94l) evaluated at go, 

AT 

= (a(g ),g Q ) = E 7o^ (a(g ), Mgo))- (96) 

v=2p+l 

However, the condition (|96|) is nugatory if the vector a(qo) is orthogonal to ker(M(g)). There 
are thus at most N — 2p parameters 70^ available for the N — 1 independent initial values 
obtainable from the initial velocity components g^, A = 1,2, . . . , N. This means that the system 
of equations fl95|) . considered as equations for the determination of the parameters jo u , v = 
2p + 1, . . . , iV is over-determined. These equations can be consistent only if the given initial 
values (Jq, A = 1, 2, . . . , iV satisfy at least (iV — 1) — (iV — 2p) = 2p—l > 1 consistency conditions, 
in addition to the constraint (1511) . This contradicts the requirement that one should be able to 
impose general initial values q^, A = 1, 2, . . . , N, for the velocities q A , A = 1,2, . . . , N, except 
for the restriction (13~4l) . which follows from the non-holonomic condition ([9]) at the initial value 
point q . 

It has thus been shown that that the initial values q A ,A = 1,2, . . . , N, which satisfy the 
condition implied by the non-holonomic constraint ([9]) evaluated at go, can not be freely chosen, 
apart from the restriction (1341) . if one demands that the d'Alembertian equations (1321) and the 
variational equations (fT4l) should have coincident solutions. 

There is still the possibility that ker(M(q)) contains only the zero vector. This is the case 
if the matrix M is regular, which is possible only if iV is an even integer. Then there is only 
the trivial solution go = 0, which naturally is in conflict with the requirement that it should be 
possible to choose the initial values g A, A = 1,2, . . . , N, freely, except for the condition implied 
by the non-holonomic constraint ([9]) at the point go- 
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5.6.2 The case ||M(g)g|| + 0. 



It has been shown above that the case (1901) necessarily leads to contradictions. It remains to 
consider the case fl9T|) . Then the conditions (l9~2l) must also be in force. 

In view of the conditions (1521) . the vector a(g) must lie in the subspace spanned by the 
basis- vectors b^q), for fi — 1, 2, . . . , 2p, 

p 

a (tf) = Z 62/*-i(g) + &2/,(g)} , (97) 

where for later convenience we have separated the sum in the representation (1971) into two parts, 
as sums over even and odd indices, respectively. The coefficients oj 2m _i and w 2ai , respectively, 
for ii = 1, . . . ,p, in the representation ( 197|) . are determined by the vector a(q). 

I will now consider such solutions z, say, of the equations (159|) which are in the linear span 
of the basis- vectors b^(q) for /i = 1, 2, . . . , 2p. Thus, 

p 

z = Z {«2 M -i&2 M -i(g) + «2 A1 &2 A1 (g)} • (98) 

H=l 



l,2,...,iV, 



Inserting g = z in Eqns. (15^1) one obtains the following equations for the components z A , A 

<q)\\ 



Y d M AB {q)z B = ^^a A {q), A = 1,2,...,N. (99) 



B=l 

At 



Contracting the equations (199|) with the basis- vector components b 2h , (g) and using the equations 
(|87|) . one obtains the following result, 

^(gX&a^iGz),*) = (M?)> a(g)), !/ = 1, . . . ,p. (100) 

Likewise, contracting the (19"5)) with the basis vector components & 2 j,_! A (g), an d using the equa- 
tions (ISol) . one finds a similar expression, 

n v {q){b 2v {q)^) = -^TTTT (^-i(g), a(g)), i/ = 1, . . . ,p. (101) 

Mini 

The equations ( llOOp and ( 11011) mean that the equations ( 159]) for the vector z have been solved 
in the sense that the equations (llOOp and (jlOip express all the components aiv-\ and a 2i ,, 
respectively, in the representation (198]) for z, in terms of one unknown scalar quantity, namely 
||M(g)z||, and the known components u 2u ~i and u 2u , respectively, in the representation (197]) 
for the vector a(g). 

The result obtained above can be summarised as follows, 

z A = ll ^ ) )| | ll E^ 1 (g) {(^(g),a(g))6 2 ^(g) - (b 2 ^(q), a(g)) & 2 /(g)} . (102) 
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It should noted that the vector z defined by Eqns. f 1 1 2 j) automatically fulfils the following 
orthogonality condition, 

(z,a(g)) = 0. (103) 

So far only the solutions z A , A = 1, 2, . . . , N, of the equations fl59|) in the linear span of the 
basis vectors b^(q), fi = 1, . . . , 2p, have been considered. The general solution q of Eqns. (E2J) 
is a sum of the solution z given above in Eqns. (I102p . and a general solution h, say, of the 
corresponding homogeneous equations. Thus, 

TV 

u=2p+l 

where the parameters 7^, v — 2p + 1, . . . , iV are free parameters. 

The general solutions of Eqns. fl59|) for the components q A are the following, 

^ : = z A + h A , A = 1,2,...,N (105) 

It should be noted that 

||M(g)g|| = \\M(q){z + h)\\ = \\M{q) z\\. (106) 
One thus finally obtains the following expressions for the general solutions of Eqns. fl59|) . 



N 

+ E 7* &/(?), A =1,2,..., TV. (107) 

i/=2p+l 

In view of the conditions ( 192]) and (I103p . one finds that the final expressions f llQ7p for the 
components q A , A = 1, 2, . . . , N, automatically satisfy the non-holonomic condition OH]). 

I now apply the results above for the initial value go at go- 

It follows from Eqns. fll07p that 

1 1 M(qo) go 1 1 p 

90 1 = I, 7~TjT E K * 1 (9o) {(M9o), a(g )) b 2l/ _ A (q ) - (b 2u -i(qo), a(q )) b 2v A (q )\ 
H a l9oJ|| u=1 

N 

+ E 70,6/(90), A = 1,2,...,N. (108) 

u=2p+l 

When the initial velocity components q A ,A= 1,2, ... ,N are given, the only unknown quantities 
in the equations f )108p are the N — 2p parameters j 0u , u = 2p+l,...,N. 

When 2 < 2p < N, there are altogether N—l independent equations among the iV equations 
in ( I108p . to be solved for the iV — 2p parameters 70^, v = 2p + 1, . . . , N, when the components 
q A , A = 1,2, ... ,N are given. This is an over-determined system of equations for the parameters 
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Jou, f — 2p+l, . . . , N, which implies that there must be at least (N— 1) — (N — 2p) = 2p — 1 > 1 
consistency conditions among the components q A , A — 1, 2, . . . , N, in addition to the constraint 
(134]) . This contradicts the requirement that it should be possible to choose the initial velocity 
<jo freely, except for the constraint 



For even N one can have 2p = N. Then the solutions for q^, A = 1,2,. . . , N are as follows, 

p 

J2 K Z\oo) {(hu(qo),a(q Q )) b 2v _ l A {q a ) - (&2„-i(gb), a(q )) K A (lo)} (109) 



. A _ \\M{q )q 
% — 



\Hqo)\\ u=1 

The N components q A , A = 1,2, ...,iV in the equations ( 11091) . are given in terms of one 
undetermined scalar quantity ||M(g ) %\\, in addition to the other known quantities in f |109[) . 
This means that there are N — 1 constraints among the components q A , A — 1, 2, . . . , N, which 
clearly contradicts the requirement that it should be possible to choose the initial velocity q 
freely, except for the constraint ([31 



In the case | \M(q ) q \ \ ^ 0, it has thus been shown that there are at least 2p— 1 consistency 
conditions on the velocity components q$. There are also N — 2p conditions of the form (192]) 
evaluated at q , namely the following conditions, 

(K(q ),a(q )) = 0, v = 2p + 1, . . . , N, (110) 

Taken together there are thus altogether at least (2p — 1) + (iV — 2p) = N — 1 conditions on the 
initial values (q^, q^), A = 1, 2, . . . , N. This is a set of stronger conditions than the consistency 
conditions involving only the velocity components q^, A = 1, 2, . . . , N. 

It has thus been shown also in the case ||M(g)o)go|| 7^ 0, that that the initial values 
c[q,A = 1,2, . . . , N, which satisfy the condition (|34l) implied by the non-holonomic constraint 
dHD evaluated at go, have to satisfy additional consistency conditions, if one demands that 
the d'Alembertian equations (1321) and the variational equations (fl4|) should have coincident 
solutions. 



5.7 Conclusion of the proof of the incompatibility theorem 

The starting point in of the beginning of the proof in Subsection 15.21 is the following set of 
equations, 

N 

Y,M AB (q)q B = Ta A (q), A=1,...,N, (111) 

B=l 

where T is a scalar parameter. The equations (11111) are necessary consequences of the assump- 
tion that the variational equations (fill) and the d'Alembertian equations (132|) have coincident 
solutions. 

It was shown in Subsection 15.31 in the case N = 3, that the equations (11 lip do not have 
solutions q A , A = 1,2,3, satisfying the required constraint if this constraint is genuinely 
non-holonomic, i.e., neither integrable as such, nor integrable by means of an integrating factor. 
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Thus the variational equations (j!4p and the d'Alembertian equations ( 132]) are not compatible 
when N = 3. 



When N > 4, the equations f illip were transformed into the following form in Subsection 



£ M AB {q)q B = ^ { ?\t a ^^ A = 1.2 JV, (112) 

where the norms 1 1 • • • 1 1 occurring above are defined by an appropriate metric tensor g. 

For N > 4 the variational equations of motion fJEj) , and the d'Alembertian equations of 
motion ( 132]) . respectively, have been considered to be initial value problems. 

The implications of Eqns. (I112p for the initial values qo and go were analysed in Subsection 
15.61 It was shown that the equations (11121) imply that the initial values q^, A = 1, 2, . . . , N, 
for the solutions of the variational equations (1141) and d'Alembertian equations (|32|) . which 
satisfy the required constraint will have to satisfy additional consistency conditions if the 
constraint OH]) is genuinely non-holonomic. 

Thus in the cases iV > 4, the d'Alembertian and variational equations of motion are in- 
compatible in the sense that in the case of a genuinely non-holonomic constraint, they can not 
have coinciding solutions with initial velocities q^, A = 1, 2, . . . , N, which satisfy the required 
constraint (Q at qo, unless these initial velocities satisfy additional consistency conditions. 

This concludes the proof of the incompatibility theorem stated in Subsection 15.11 



6 Discussion 



In this paper I have refined and completed an earlier proof [3J of the incompatibility of the 
variational equations of motion ( 1T4"P and the d'Alembertian equations of motion (13"2"p . It was 
observed in Ref. [5], that the proof in Ref. [3] is not valid in the cases p — 1 and N > 4, 
since the conditions of consistency derived in the proof in Ref. [3] are nugatory when p = 1 
and N > 4. This circumstance has now been taken into account in the present proof, in which 
more precise consistency conditions have been derived than in the earlier proof. 

In the present proof the formulation respects explicitly the desirable circumstance that the 
geometry of configuration space need not be Euclidean, but rather e.g. Riemannian. The tensor 
machinery necessary in this case, which was developed in Subsection 15. 4\ and in particular in 
Subsection 15. 5[ may seem a little heavy. It is a price one has to pay for the desire to handle a 
Riemannian configuration space using a classical tensor formulation. 

The proof in Ref. [3] avoided using a full-fledged tensor formalism, by using essentially only 
local analysis related to the fixed initial point qo in configuration space for the cases N > 4, and 
by implicitly using so-called normal- or geodesic co-ordinates [UJ at that point. This means 
that one was permitted to use a Cartesian metric at the fixed point qo in Ref. [3]. 
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The incompatibility theorem proved in this paper does not mean that one is always un- 
able to associate some variational principle with systems with non-holonomic constraints. It 
has only been proved that the particular variational equations (|14p are incompatible with the 
d'Alembertian equations of motion when N > 4, in the sense that these equations do not have 
solutions with general initial velocities, if the homogeneous constraint used in this paper is 
genuinely non-holonomic. For iV = 3 a stronger result was proved, namely that the equations 
of motion in question do not have coinciding solutions at all, if the constraint is genuinely 
non-holonomic. 

Conversely, it is by no means guaranteed that the variational equations ( !T4|) . and the 
d'Alembertian equations (|32|) have coinciding solutions, if the initial values q^, A = 1, 2, . . . , N, 
satisfy those partly explicit and partly implicit consistency conditions which have been derived 
in this paper. For instance, the consistency consistency conditions involving only the initial 
velocities do by no means exhaust all the implications of the equations flll2p for the initial val- 
ues qo and qo, not to mention additional consistency requirements which come into play when 
one takes into account the full information residing in the equations of motion. An example 
of reasonably simple yet intricate additional consistency conditions on the initial values go was 
given for an important class of solutions of the equations (j!12p in Subsection 15.6.21 
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